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Herschel images in six photometric bands show the thermal emission of the debris disk surrounding /? Pic. In the three PACS bands 
at 70 /im, lOOyum and 160yum and in the 250 //m SPIRE band, the disk is well-resolved, and additional photometry is available in the 
SPIRE bands at 350/jm and 500 yum, where the disk is only marginally resolved. The SPIRE maps reveal a blob to the southwest of/3 
Pic, coinciding with submillimetre detection of excess emission in the disk. We investigated the nature of this blob. Our comparison 
of the colours, spectral energy distribution and size of the blob, the disk and the background sources shows that the blob is most likely 
a background source with a redshift between z = 1.0 and z= 1.6. 
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1. Introduction 

The circumstellar disk around the A6V starts Pic was discovered 



by IRAS in 1983 ( Aumann et al. |1984i and especially the in 



ner part has been studied extensively ever since (Kalas & Jewitt 
T995||Paiitin et al.J! 1997, Mouillet et al.|p97^ Heap et al.^ 2000, 



Telesco et al. 2005[ l. Optical observations imply that the disk ex- 
tends to 95" offset, which corresponds to 1800 AU ( |Larwood^ 
[Kalas |2001[ ). Warps and asymmetries in the inner part of the 
disk suggest the presence of planetesimals and even planets, the 
latter of which has been confirmed by |Lagrange et al.| ( |2010[ ). 
Owing to its relative proximity to the Earth (19.44 + 0.05 pc; 
2007 |l, the debris disk was resolved spatially at 
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longer wavelengths as well (Holland et al.[ |1998t [Nils son et al.[ 
|20TOi|Liseau"etam2003] l. 

The p Pic debris disk has been imaged by Herschel in six 
photometric bands. In the 250-500 yum bands a blob appears to 
the southwest of the star This blob coincides with the one ob- 



served at 850yum with S CUBA by Holland et al.| ( 
870 yum with LABOCA by |Nilsson et al.| ( |2009| l and 
of authors assumed it to be a structure in the disk. 
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A first counter-argument was already mentioned in [Dent 



et al. d2000jl. Assuming the blob indeed originates in the debris 



disk and has the same dust characteristics, its observed emission 
at 850jum can only be explained by at least 1.2 Mmooh, a mass 
comparable to the total mass of the disk in small dust particles. 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 
** Postdoctoral Fellow of the Fund for Scientific Research, Flanders 



Given the blob's distance from the central star and the fact that 
there is no counterpart for such a large amount of dust in deep 
optical observations, [Dent et al. (2000) argued that the blob can- 
not be a structure in the disk. Vandenbu ssche et al., ( ,2010} argued 
that given the high amount of background galaxies in the neigh- 
bourhood of y6 Pic, the blob probably is a background source. 

In this paper, we investigate the hypothesis that the blob is a 
background source and not a structure in the circumstellar disk. 



2. Observations and data reduction 

Photometric maps were obtained with the Herschel Space 
Observatory ( [Pilbratt et aTj [20T0) l with the PACS and SPIRE 
photometers. A concise observing log can be found in Table [T] 
More information concerning the scientific capabilities, obser- 
ving modes, data reduction, calibration, and performance of 
these instruments can be found in |PogIitsch et al.| (|20I0 1 for 
PACS, and in [Griffin et"aL] ( [20T0| ) and [Swinyard et al.| ( |20l"0 1 
for SPIRE. 



The observations were performed in the context of the 
'Stellar Disk Evolution' Guaranteed Time Key Programme (PI: 
G. Olofsson) during the science demonstration phase. 

Pipeline processing as well as image analysis have entirely 
been made in the interactive processing environment HIPE, 



which has been developed for the Herschel project (Ott & 
[Science Ground Segment Consortium! |2010| l. 
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Figure 1. Surface brightness maps in the SPIRE bands (in the upper panels) with the fitted Gaussian overplotted; and the maps 
after subtraction of the fitted disk, with the position of the blob indicated with a 20" radius circle in the lower panels (left: 250 /zm; 
middle: SSO/zm; right: 500 yum). All maps display the same area on the sky and are centred on the optical position of /? Pic. 



Table 1. Observation log for the photometric Herschel obser- 
vations of yS Pic. Column 1 gives the instrument, column 2 the 
observation identification, column 3 the date of the observation, 
column 4 its duration and column 5 the observed wavelengths. 





Observation 


Date 


Duration 


Filters []im\ 


SPIRE 


1342187327 


2009-11-30 


3336 s 


250, 350, 500 


PACS 


1342185457 


2009-10-07 


808 s 


100, 160 


PACS 


1342186612 


2009-11-01 


5506 s 


70, 160 


PACS 


1342186613 


2009-11-01 


5506 s 


70, 160 



2.1. SPIRE observations 

The SPIRE observation includes three wavelength bands simul- 
taneously (at 250 /vm, 350 /vm and 500 /im), covering an area of 
8' X 8' homogeneously. The maps were constructed with the 
naiveMapper task. This algorithm projects each sample onto a 
single map pixel. The pixel size is 4", 6" and 9" in the three 
bands respectively, and the beam FWHM is 18.2", 24.9" and 
36.3" respectively. 



2.2. PACS observations 

At 100 jum, we have a single shallow PACS scan map, taken at 
medium speed (20"s"'); while at 7Qfj.m and 160 /^m, a scan and 
a cross-scan at low speed (10"s"') are at our disposal. For the 
shallow observation, an area of 2'x2' is covered homogeneously; 
for the deep observations an area of 2.5' x 2.5'. 



The PACS observations were processed with the HIPE 
pipeline script for the reduction of observations of extended 
sources, such as /3 Pic. Removal of cosmic rays was performed 
with the IlndLevelDeglitch task. For each pixel, cr-clipping was 
applied on all contributions before projecting the remaining val- 
ues onto a map. The flux conversion was made with version 6 
of the response calibration. Before applying a high-pass filter to 
eliminate low-frequency drifts, the region around the source was 
masked. The width of the filter window was taken as half the 
scan leg length (3.7"). We used the photProject task to combine 
the detector time series signals into maps of 1" pixels in the blue 
and green, and 2" pixels in the red, weighting the signal by the 
respective noise. 

3. Image analysis 

In this section we discuss how we determined the flux of the 
debris disk, the blob, and the background sources. 

From the SPIRE maps shown in Fig.fT] a blob to the south- 
west of the debris disk of y6 Pic can be clearly seen (at an offset 
(Aa, A5) = (-25", -23") from the star), after subtraction of the 
flux of the disk. In Fig.l2] we notice tens of background sources 
in the 10' x 10' region around the star We quantified whether 
the far-IR colours of the blob resemble the colours of these other 
sources more than the colours of the debris disk itself. If so, this 
would be a strong indication that the blob is a background source 
and not a structure in the disk. 

In the PACS maps, the background sources are less promi- 
nent, and the complex shape of the PACS point-spread-function 
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Table 3. Flux of the debris disk for PACS and SPIRE (calculated 
for a 60" target radius). 



Figure 2. SPIRE map at 250 //m on which the background 
sources are indicated with green circles. The green rectangle in- 
dicates the region that is covered by the two scan directions. 

Table 2. Fit parameters for the subtracted Gaussian disk model 
in the SPIRE wavelength bands. 



Parameter 


250 yum 


350/im 


500 //m 


a,,„„r [h m s] 


05 47 17.12 


05 47 17.11 


05 47 16.94 


Ocenler V J 


-5103 55.5 


-5103 55.52 


-5103 58.20 


FWHM™,„, ,„,["] 


23.20 


28.96 


43.01 


FWHM„,„„ „,,,["] 


17.76 


23.54 


34.78 


Position angle [°] 


37.40 


39.31 


39.31 



(PSF) prohibits a swift separation of the disk flux from the PACS 
maps. Moreover, the surface brightness profiles for PACS in 
Vandenbussche et al. ( 2010, see Fig. 3) suggest the blob is faint 
and cannot be discerned in the PACS wavelength range. 

In the SPIRE maps, we separated the fluxes of the blob from 
the fluxes of the debris disk by fitting an elongated Gaussian to 
the observed disk and subtracting it, using the sourceFitting and 
imageSubtract tasks in HIPE. Fig.fT] shows the original SPIRE 
maps on which the respective fitted Gaussians were overplotted, 
and the residual maps with the blob indicated. The fit parameters 
are summarised in Table |2] 

We can now determine the flux of the debris disk and of 
the blob by performing aperture photometry on the original 
maps and on the difference maps with the rectangularSkyAper- 
turePhotometry task in HIPE. This task calculates the flux within 
a circular aperture and estimates a sky contribution from a rect- 
angular region close enough to the target to achieve a similar 
coverage. 

To calculate the flux of the disk, we used a circular 60" ra- 
dius target aperture. A value for the sky contribution was ob- 
tained with a sky estimation algorithm, corresponding to the 
mmm routine from the IDL Astronomy Library to estimate the 
background in a stellar contaminated field. The acquired fluxes 
are summarised in Table[3] The errors on the fluxes are do- 
minated by uncertainties in flux calibration, which are less than 
15% for SPIRE ( ,Swinyardetal.|[20T0: ), less than 10 % at 70 /zm 
andlOO/m, and less flian 20% at 160//m for PACS ( |Poghtsch 
|et al.[|20T0| . To obtain the flux of the blob, we used a target ra- 
dius of 20" on the difference maps. The results are summarised 
in Table |4] In the difference maps, the residuals from the disk 
subtraction become substantial, especially at 250 /vm. The errors 
on the fluxes of the blob are hence calculated as the standard de- 
viation on the flux in 20 circular apertures (with a radius of 20") 
close to the subtracted disk. 

To determine an upper limit for the flux of the blob in the 
PACS maps, we calculated the flux within a 10" aperture at the 



Wavelength [}im\ 


Flux [Jy] 


70 


15.34 ±0.8 


100 


9.34 ± 0.5 


160 


4.55 ± 0.5 


250 


1.88 ±0.14 


350 


0.81 ±0.06 


500 


0.37 ± 0.03 



Table 4. Upper limits for the flux of the blob in the PACS images 
and the flux of the blob in the three SPIRE bands (calculated for 
a 20" target radius), and the corresponding beam FWHM. 



Wavelength [jum] 


Flux [Jy] 


Beam ["]" 


70 


<0.39 


5.26x5.61 


100 


<0.16 


6.69 X 6.89 


160 


<0.15 


10.46 X 12.06 


250 


0.081 ±0.027 


18.2 


350 


0.050 ± 0.006 


24.9 


500 


0.026 ± 0.004 


36.3 



Notes. <"> Taken from the PACS Observer's Manual and SPIRE 
Observers' Manual (version 2.4). 



position of the blob (in the original images) and subtracted the 
flux within a 10" aperture on the opposite side of the star Scaling 
the difference in flux between the SW- and NE-aperture with the 
encircled energy function yields the upper limits listed in Table 

a 

To identify the background sources, we applied the HIPE 
sourceExtractorSussextractor task to the original map at 250 jum 
on the homogeneously covered region (with signal-to-noise 
threshold of 5.0). This task yields a list of sky coordinates and 
fluxes for the extracted sources, which can be fed into the task 
as priors when applying it to the two other original maps. This 
forces the task to extract fluxes for the same sources in all three 
images. Sources within a radius of 50" from the star were re- 
moved from the list and 50 sources were left for analysis (see 
Fig.l2|. Their sky positions and fluxes in the SPIRE bands are 
summarised in Table |6] 



4. Discussion 

4. 1 . Comparison of the far-IR colours of the blob, the disk 
and the background sources 

Fig.|4]shows /^5oo/^35o vs. F250/F3SQ colour-colour diagram for 
the debris disk, the blob and the background sources. The 
colours of the blob are more similar to the colours of the back- 
ground sources than the colours of the disk. This is a first indi- 
cation that the blob is a background source and not a structure 
in the disk, but no decisive conclusion can be drawn due to the 
relatively large error bars. 

4.2. Comparison of the size of the blob and the background 
sources 

A second argument that the blob is a background source was 
found in the comparison between the size of the blob and the 
size of the background sources. The full-width at half-maximum 
(FWHM) of the background sources was determined in the three 
SPIRE bands by fitting a Gaussian to each of the sources in the 
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7Ufim 

0,00 




Figures. Surface brightness maps in the PACS bands with the 10" 
blob overplotted (left : lOfim; middle : IQQfim; right : 160 yum). 



apertures used to determine an upper limit for the flux of the 
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Figure 4. Colour-colour diagram of the /3 Pic debris disk, the 
blob and the background sources. The error bars for the fluxes 
of the background galaxies are not visible on this scale. 

Table 5. Median FWHM of the background sources and the 
FWHM of the blob in the three SPIRE bands. 



Wavelength [jum] 


Background sources ["] 


Blob ["] 


250 
350 
500 


18.7 ±5.5 
25.7 ± 11.2 
38.2 ± 19.5 


23.4 
33.4 
37.2 



original maps. A Gaussian fit to the blob in the difference maps 
yielded the FWHM of the blob in the three bands. The results 
are listed in Table|5] 

The blob is not resolved, and its size is clearly consistent 
with the size of the background sources. This is a second indica- 
tion that the blob is a background source. 



with fluxes obtained with Herschel, SCUBA, LABOCA, and 
SIMBA, normalised to the flux of the disk at 250 yum, as well as 
the flux points for the blob, normalised to the flux of the blob at 
250jum (in red). Overplotted are three modified Rayleigh-Jeans 
laws {Fy ex i/+^). The spectral index fi describes the mean dust 
opacity {k cx v^) and is an indicator of the grain size distribution 
of the dust. The SED of the disk is best described by a spec- 
tral index/? - 0.34 (implying (sub)yum grains), while a negative 
spectral index is needed to fit the SED of the blob (/? = -0.45 
produces the best result). The latter can only be explained by 
a steep temperature gradient that cannot occur on the scale the 
blob would have if it were a structure in the disk (i.e., about 
470 AU in diameter). 

The submillimetre fluxes of the blob and the disk, as detected 
by SCUBA, LABOCA, and SIMBA, are summarised in Table|7] 

4.4. Determination of the redshift 

Assuming the blob really is a background source, we can es- 
timate its redshift based on the far-infrared fluxes. Amblard] 
et al. ( 2010| l investigated the redshift distribution of submil- 
limetre galaxies in the Herschel- ATLAS survey. Their ^500/^350 
vs. F25q/Ft,sq diagram shows the average redshift in that colour- 
colour space. In Fig.l6] we have overplotted the datapoint for the 
blob. Within the error bars, we find redshifts between 0.78 and 
2.55. 

Assuming the blob is a background galaxy with redshift 
z - 1.6, we derive a luminosity of the order of 10^^ W Hz ' at 
250//m (from the integrated flux at that wavelength). The lumi- 
nosity distance di we used for this calculation was derived from 
the relationship in Fig. 1 o f |Choudhury & Padmanabhan ( 2005| l. 
The luminosity function of |Eales et al. '{ 2010 their Fig. 4) shows 
this is an unlikely high value. Higher redshifts would result in 
even higher luminosities. Because this is improbable, we esti- 
mate z to be no higher than 1.6. 

Rejecting the 5% lowest values for the redshift within the 
overplotted box (see Fig. IT] for the histogram), we estimate z to 
be at least 1.0. 



4.3. Far-infrared SED of the blob and the disk, and 
interpretation as dust spectrum 

Another test for the hypothesis that the blob is a background 
source are the relative spectral energy distributions (SEDs) of 
the disk and the blob. Fig. |5] shows the SED of the debris disk 



5. Conclusions 

We have investigated the nature of the blob that has been de- 
tected to the southwest of the debris disk of/? Pic not only with 
Herschel, but also with SCUBA, LABOCA, and SIMBA. By 
comparing the submillimetre colours of the blob with the colours 
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Figures. Spectral energy distribution of the debris disk (in 
black) and the blob (in red), normalised to their repective flux 
at 250 jum. The Herschel fluxes of the disk were obtained by 
integrating over a target radius of 60", the SPIRE fluxes of 
the blob by integrating over a target radius of 20". The PACS 
fluxes of the blob are upper limits, determined as described 
in Sect.[3| Add itional flux points for the disk were taken from 
HoUand et al!] ^^ (S CUBA: 850 / /m), Nilsson et al] ^Jm\ 
(LABOCA: 870 //m) and |Liseau et"aL|(2003) (SIMBA: 1.2 mm). 
The SCUBA and SIMBA fluxes were integrated over a 40" ra- 
dius aperture while the LABOCA flux point was obtained by 
fitting a Gaussian to the disk. Additional flux points for the 
blo b were taken from Dent et al. ("2000) (SCUBA : 850yL(m) 
and [Nilsson eFaT] (12009) (LABOCA : 870 yum). Overplotted are 
three modified Rayleigh-Jeans laws {F^ oc i/+^). The spectral 
indices have been determined using the fluxes at wavelengths 
A. >160/im for the disk and A >250/vm for the blob, with the 
inverse of the errors on the fluxes used as weight. Note the large 
error bar for the 870 /im flux point of the blob. The /? = curve 
is shown for comparison. 



of the debris disk itself and with the background sources that 
are present in the images, we conclude that the blob is not a 
structure in the disk, as previously suggested, but a background 
galaxy. Support for this hypothesis was found from differences 
in the spectral indices and SEDs between the disk and the blob, 
and the fact that the size of the blob is consistent with the size 
of the background sources. Comparing the colour of the blob to 
the colour-colour diagram from |Amblard et aT] ( |2010| l yielded a 
rough estimate for the redshift of the blob {z between 0.78 and 
2.55), which we narrowed down to z between 1.0 and 1.6. 

Ultimately, the nature of the blob could be probed by spec- 
troscopy. The redshift of a galaxy can be determined based on 
the CO emission, because the spacing between two successive 
transitions shrinks with redshift (as z + 1). The large total band- 
width and high angular resolution make ALMA the ideal instru- 
ment to study the CO emission of the blob, and determine its 
redshift. For a normal galaxy like the Milky Way, with z = 3, 
the study of the CO emission takes less than 24 h observing time 
( IWootten & Wilson|[2005| ). 

Band 3, the most sensitive available band, will show CO 2-1 
transitions for z - 1-2. Higher (lower) redshifts would move 
CO 2-1 out of the band, and CO 3-2 (CO 1-0) into the bandpass. 
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Figure 6. Colour-colour diagram of the Herschel-P^l^AS 
sources ( Amblard et al.[ 2010) 1. The coloured background in- 
dicates the average redshift for 10^ randomly generated model 
SEDs. The SPIRE colours of background sources can be over- 
plotted to determine their redshift. Other object types (e.g. the 
debris disk itself) should not be overplotted. The small crosses 
are the 1947 galaxies detected in the context of the ATLAS 
open-time key programme (with a significance greater than 5cr 
at 350 yum and 3cr at 250 and 500 jum). The overplotted rectangle 
indicates the position of the blob in this colour-colour space. The 
redshift in this rectangle lies between 0.78 and 2.6. 
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Figure 7. Histogram of the redshift distribution in the rectan- 
gle marked in Fig.|6] Values range between 0.78 and 2.55. The 
shaded region indicates the redshift to which we have narrowed 
down the region. 



Observing additional lines in, e.g., band 7 (with the penalty of 
lower sensitivity), or JVLA (for CO 1-0) can rule out any ambi- 
guity. 
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Table 7. Submillimetre fluxes of the disk and the blob as de- 
tected by SCUBA, LABOCA, and SIMBA. The blob was not 
observed at 1 .2nim 



wavelength [}im\ 


flux disk [mjy] 


flux blob [mJy] 


References 


850 
870 
1200 


104.3 ± 10.0 
63.6 ± 6.7 
35.9 + 9.7 


24 ±5 
7.0 ± 14.9 

/ 


1,2 
3 
4 



References. (l) |Holland et"aL](T998l >; (2) |Dent et al.]pOOOl ); (3) |Nilsson| 
[etaLlp009t ; (4) |Liseau et aL|p003) 
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Photometry of the background galaxies in the SPIRE 
wavelength bands 

Table [6] summarises the positions and photometry in the three 
SPIRE bands for the extracted background galaxies. 

Submillimetre fluxes of the disk and blob 

In Table |7] the submillimetre fluxes of the disk and the blob as 
measured wifli SCUBA, LABOCA, and SIMBA are listed. 



S. Regibo et al.: A background galaxy in the field of the/? Pic debris disk (RN) 
Table 6. Positions and photometry of the background sources in the SPIRE wavelength bands. 
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